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ABSTRACT 

Many  times  populations  exist  which  logically  must 
satisfy  a  stochastic  ordering  requirement.  Nevertheless, 
estimates  of  these  populations  may  not  bear  out  this 
stochastic  ordering  because  of  the  inherent  variability  of 
the  observations.  This  paper  will  consider  the  problem 
of  finding  maximum  likelihood  estimates  of  stochastically 
ordered  survival  functions  for  the  cases  a)  one  survival 
function  being  fixed  in  advance  and  b)  estimating  both 
survival  functions  when  the  data  includes  censored 
observations . 

A  numerical  example  is  handled  in  detail  to  illustrate 
the  solution  to  this  problem. 

AMS  Classification  numbers:  Primary  60G05;  Secondary  62N05. 


Key  Words  and  Phrases:  maximum  likelihood  estimation; 
survival  functions;  stochastic  ordering;  censored 
observations;  Kaplan-Meier  product  limit  estimator; 
order  restrictions. 


I.  INTRODUCTION 


Many  times  a  new  population  is  created  which  can 
logically  be  only  stochastically  greater  (or  less)  than 
the  old  population.  Nevertheless,  estimates  of  these 
populations  may  not  bear  out  this  stochastic  ordering 
because  of  the  inherent  variability  of  the  observations. 
Brunk,  et. al.  (1966)  have  given  Maximum  Likelihood 
Estimates  (M.L.E.'s)  of  the  C.D.F.'s  of  two  stochastically 
ordered  distributions  when  all  observations  are  complete 
and  the  estimated  distributions  are  required  to  be  of  the 
discrete  type.  This  paper  will  consider  the  similar 
problem  of  finding  M.L.E.'s  of  stochastically  ordered 
survival  functions  for  the  cases  a)  one  survival  function 
being  fixed  in  advance  and  b)  estimating  both  survival 
functions  when  the  data  includes  censored  observations. 

The  unrestricted  (except  for  the  requirement  of 
being  a  discrete  distribution)  M.L.E.  of  the  survival 
function  when  censored  data  is  involved  has  been  developed 
by  Kaplan  and  Meier  (1958).  Since  the  results  of  Kaplan 
and  Meier  are  so  well-known  and  widely  used,  this  paper 
will  largely  conform  to  the  notation  developed  there. 
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II.  NOTATION,  THE  PROBLEM,  AND  THE  SOLUTION 

Initially,  we  consider  the  following  problem.  Inde¬ 
pendent  observations  are  taken  from  a  discrete  distribu¬ 
tion  on  the  positive  part  of  the  real  line  with  survival 
function  P(t).  We  wish  to  find  the  M.L.E.  of  P(t)  subject 
to  P(t)  s  Q(t)  for  all  t  where  Q  denotes  the  survival 
function  of  a  fixed  discrete  distribution  with  only  a 
finite  set  of  points  possessing  positive  probability. 

Complete  observations  (which  we  will  call  deaths)  occur  on  a 


subset  of  the  times  S.  <  S~  <  •••  <  S  (Sn  =  0  and  S  ..  = 

i  i.  mu  m+l 

convenience) .  The  number  of  deaths  at  Sj  is  6j .  We  let 
Aj  denote  the  number  of  censored  observations  (losses) 
in  [  Sj  ,  Sj  +  1),  assumed  to  occur  at  L^^,  i  *  I,***,  • 

In  Section  III,  we  will  prove  that  £(t)  may  be  expressed 


in  the  following  manner.  Let  <  S2  <  . . .  <  Sm 
denote  the  ordered  values  of  the  times  of  death  combined 


with  the  points  of  positive  probability  under  Q,  and  let 
m 

ni  *  j=^  Jj  *  *j  denote  the  number  of  items  surviving  just 

prior  to  .  Our  Pj  and  qj  are  related  to  P(*)  and  Q(*) 
respectively  by 

Pj  =  In  [P(Sj)/P(Sj_1)]f  and 
q.  =  In  [(KSjj/QtS.^)]. 

Then  the  restricted  M.L.E.  for  t  <  S  is  given  by 

m 

P(t)  -  exp [  Z  p.] 
i;Sist 

where  p^  is  given  in  the  following  theorem. 

(Throughout  the  paper  we  shall  treat  00  and  -°°  as  real 
numbers  greater  than  and  less  than  all  other  real  numbers 
respectively.  We  shall  also  adopt  the  conventions  that 
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fcn(O)  =  ,  +»/+<»  =  +1  ,  0 (+*>)  =  0  ,  O/O  *  1  ,  and  0  = 


Theorem  2.1.  Let  k + .denote  the  constant  k  such  that 
• -  a  ,b 


(2.1) 


b  n  •  -  5 .  +  k  b 

f  <nt  3  «: ■  ?  ij- 


a  -j  a 

if  k  >  0  and  exists  and  0  otherwise.  Then  if 

min  max  . 


(2.2) 

A 


A 

^i  -  asi  i<b  ka,b  * 


p^  is  expressable  as 


n-  -  5-  ♦  k- 

A  n  f  1  1  1 

p_  =  fn  ( - 7T - ) 


ni  +  ki 


Of  course  if  the  last  observation  at  say  t*  corres¬ 
ponds  to  a  loss,  the  P(t)  may  be  defined  arbitrarily  for 

A  A 

t  >  t*  providing  P(t)  2.  Q(t)  and  P(t)  is  a  survival  function. 

(Note  that  if  Q(S1)  -  Q(S2)  -  •••  ■  Q(Si)  =  1,  then 

i  n .  -  S .  +  k  i 
l  tn(  )  -  Iq,  *  0 

has  the  solution  k  =  °°.  This  leads  to  the  intuitively  reason 
able  solution  p^  *  p2  =  ***  =  p^  =  0,  even  though  the  likeli¬ 
hood  is  identically  zero  in  this  case.) 

Alternatively,  the  p ^ ' s  can  be  found  more  easily  by  the 
following  algorithm: 

1.  Find  the  largest  k^  >  0  such  that 


n.  <5.  +  k 


1  tn  t  3n  I  V  '  )  '  1  qj 
1  ^  1  1  J 


for  0  <  i^  < 


If  more  than  one  i^  works,  choose  largest. 

a  n.  6.  +  k. 

2.  Let  P.  ■  tn  (-3— ^ — *-r — -)  for  1  s  j  s  i 

3  nj  +  1 


3.  Find  the  largest  k2  >  0  such  that 

i2  n,  6.  k,  12 

I  In  (-J— ^ ^ r — -)  =  £  q.  for  i,  <  i,  «  m. 

i1+l  nj  k2  ix+l  3  1 

If  more  than  one  i2  works,  choose  largest. 

n .  -  6 .  +  k2 

4.  Let  Pj  =  In  (  Jn.  - ),  ix  <  J  "  l2 ' 

J  L 

5.  1-tc.  If  at  some  point,  no  such  positive  k( 
exists,  then 

n .  -  6  . 

p.  =  tn(-J - i-)  for  i^_1  <  j  m. 

3  n . 

J 

/V 

If  the  last  observation  at  say  t*  is  a  loss,  then  P(t)  may  be 

may  be  defined  arbitrarily  beyond  t*  as  long  as  P (t)  "  ^(.ti 

/\ 

and  P(t)  is  a  surivial  function. 

The  order  constraints  P(t)  s  Q(t)  can  be  handled  similarly 

with  the  modifications  given  below. 

If  6.  =  6  =  *  *  '  =6 .  =  0  and  6.  >  0,  set  p.  =q. 

12  l  l  +1  ’  J  J 

o  o 

for  j  =  l,2,‘**,i  .  I'or  j  >  iQ  ,  Theorem  2.1  will  still 

work  by 

i)  restricting  a  to  be  greater  than  i  , 

ii)  defining  ka  ^  to  be  the  value  k  in  (2.1)  if  k  -  0 
and  exists,  and  0  otherwise,  and 

iii)  interchanging  the  words  min  and  max  in  (2.2). 

The  algorithm  is  modified  to  handle  P(t)  s  Q(t)  by 


i)  setting  Pj  =  for  j  -  1,  2,  **’,  iQ  , 

ii)  beginning  the  sum  in  step  1  at  i  +1  rather  than  1, 

iii)  replacing  k .  .■>  0  by  k.  <  0  ,  and 


s 


ivj  replacing  "largest  >  0"  by  "smallest  k^  0"  in 

steps  1  and  3  and  replacing  "positive"  by  "negative" 
in  step  5. 

The  two-sample  problem  where  a  second  set  of  independent 
observations  are  taken  from  a  distribution  with  survival  func¬ 
tion  Q(t)  can  be  handled  by  essentially  the  same  methods  already 

given.  In  this  case,  let  S.<  ’**<  S  denote  the  ordered  times 

1  m 

of  death  of  the  combined  samples.  We  let  m^  denote  the 

number  of  items  in  the  second  sample  surviving  just  prior 

to  S.,  and  d.  the  number  of  deaths  at  S.  in  the  second'  sample. 

J  J  J  ' 

Then  Theorem  2.1  will  still  apply  if  (2.1)  is  replaced  by 

k  n.  -  6.  +  k  b  m.-d.-k 

U.i  )  t  fa(-L-  Jy  )  .  j  tn  (-L-  j  y  ) 

a  j  a  ,,lj  * 

subject  to  a  special  case. 

As  is  to  be  expected,  when  Theorem  2.1  applies 

A 

.  m .  -  d .  -  k . 

-  «»(-> - V-1)  • 


■j 


m  .  -  ft  . 

J  i 


[f  d.  =  d0  =  •••  =  d.  =  0,  d.  x1  >0,  then  we  must 
12  2o  lo 1 

set  Pj  =  ,  j  <_  ic,  in  our  likelihood  function.  Tt  is  straight 


forward  to  show  that. 


(2.31 


m  .  +  n  . 

-J _ 1 

m.  +  n 


are  the  solutions  for  p^  and  q^  in  our  equations.  The  earlier 
scheme  then  works  if  we  require  a  to  exceed  i 

0  ’ 

The  algorithm  will  still  work  if  we 

i)  define  and  q^  as  in  (2.3)  if  j  <  iR 
ii)  replace 


•  E  qi 


by 


E 

^-1+1 


m. 

£n(-^- 


d. 

—L 


-)  for  all  t 


m. 

1 
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iii)  and  define 

m.  -  d-  -  k. 

1j  ‘  *n(  m.  -  k, - >  f°r  H-l  <  3  i  %  • 

3  *■ 

III.  DERIVATION  FOR  THE  ONE  SAMPLE  PROBLEM 


Clearly,  the  likelihood  function  of  the  observations 
is  expressible  as 


o  (0)  m  6.  Xj  m 

(3.1)  L(P  (t) )  =  II  P(L.  )  •  II  {(P(S.  -  0)-P(S.)]  3  n  P(L,PJ)} 

i=l  1  j=l  3  3  i=l  1 

where  S.  <  S_  ...  <  S  are  defined  as  in  Section  II.  We  wish 
li  m 

to  maximize  this  expression  subject  to  the  condition  that 
P(t)  >  Q(t)  for  all  t.  Clearly  to  do  this,  the  P(L^-^)  and 
P(Sj  -  0)  should  be  as  large  as  possible  and  the  P(Sj)  as  small 
as  possible.  Decreasing  the  P(Sj)  too  much,  however,  may  cause 
the  order  constraints  to  be  violated.  However,  since  we  are 
dealing  with  discrete  distribution,  clearly  it  will  suffice  to 
maximize 


m  6,  A. 

n (P(S.  ,)  -  P(S.)]  3  PCS.)  3 

J  J  A  J  J 

subject  to  P(Sj)  >  Q (S j ) ,  j  =  0,  ...,  m,  since  we  can 

take  P  (■* )  to  be  constant  on  the  intervals  [S.,  S^  +  ^). 
Equivalently,  we  wish  to  maximize 


m 

n  [l  - 

l 


PCS.)  6,  PCS.  .) 


pTs 


j-D 


j  '  2J 


6i  PCS  ) 

'FTSjtj)  '"P(V 


PCs.)  Q(s . ) 

or,  letting  p'  =  and  qj  =  , 

to  maximize 


m  ft ■  A .  A . +6  . 

n  1(1  -  P-)  3  3  ’  n  pT  ^  •’l 

j  =  l  3  3  i<j  1 


] 


A  . 

J  » 


subject  to 

i  i 

II  p'  i  II  qT  for  i 
j=l  3  j=l  3 


1,2, 


.  . . ,  m  and  1 


0  for  all  i  . 
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m 

If  we  let  *  E  5.  +  A.  denote  the  number  of  items 
j  =  i  J  J 

surviving  just  prior  to  Si,  our  expression  becomes 
m  6  .  n . -6  . 

H  (1  -  Pr)  J. 

Finally,  making  the  change  of  variables, 
pi  =  In  p^  and  q.^  =  In  q' 

and  considering  the  natural  log  of  the  likelihood,  our 
problem  is  to  maximize 

m 

(3.2)  f  (p  ,  .  .  .,p  )  =16.  fn(l 
1  "l  l  J 

subject  to  the  constraints 

i  i 

I  -■  j  q.  and  0  i  p.  >  -  *  f0r  i  =  1(  2,  .  .  .  ,  m. 


Suppose  the  constraint 
i  i 

Ep.  =  Eq.  =  c(0  >  c  >  -  “>) 

1  J  1  J 

is  imposed  when  maximizing  f(pj,  ...»  Pm)  .  Then  by 

i-1 

writing  p.  =  c  E  and  setting  the  partial  derivatives 

1  l  Pj 

equal  to  zero,  we  obtain  the  system  of  equations 
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Solving  these  equations  gives  the  values 
n .  -  6  . 

Pi  *  fn(-2_ - ■*-) »  J  >  i»  and 

J  j 

n .  -  6 .  +  k 

Pj  ■  <n(  Jnj  *  Ic - >•  i  K 

(Strictly  speaking , cquat i ons  (3.2)  ami  (3.5)  are  not  valid  ii 
6j  =  0.  However,  our  solutions  are  still  correct 
according  to  our  previous  conventions.)  Moreover,  noting 
that 

«  ,ni  ■  si  *  k. 

Pi  ■  ln(  F-.-E - >* 

by  adding  the  first  i  equations  it  easily  follows  that  k 
must  be  a  real  solution  of  the  equation 

i  n.  -  6 .  +  k  i 

1  — rh? — )  *  s  <ii  *  c. 

i  nj  K  i  j 

In  most  situations,  k  will  not  have  a  closed  form 
expression. 

If  m*  constraints  are  imposed  in  obtaining  the  M.L.E., 

say 


l  * 

m* 


xm* 


*1  11 

Z  p.  =  q  =  ci>  •••»  £  P .  ^  H  4  cm*» 

1  J  1  J  1  1  j  1  J  m 


then  this  is  equivalent  to  the  constraints 


m' 


^  Pi  =  ci»  ^  Pi  =  c2  cl*  •••>  .  *•  Pi  ~  cm*  '  cm*. 

1  J  1  i.+l  3  L  1  i  *+l  J  m  m 

1  m* 
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Each  part  of  f ( * )  may  then  be  maximized  separately 
to  obtain  a  solution  of  the  form 
A  n.  -  6.  +  k„ 

(3.4)  Pj  =  - ).  if  <  j  *  V  l  -  1,  ....  m*  +  1. 

where  k^  is  the  unique  real  solution  of  the  equation 

(3.5)  it  n,  -  6,  +  k, 


£n( 


h-i*1 


j - i_ 

n!  +  k  / 


-) 


°t  '  Cl-\*  ^  "  i >  m** 


(We  take  iQ,  cQ,  and  km*+1  to  all  be  zero.) 

The  problem  of  finding  the  stochastically  ordered 
M.L.E.  of  the  survival  function  is  thus  that  of  determin¬ 
ing  which  constraints  should  be  imposed. 

To  answer  that  question,  the  following  lemmas  are 
important . 

Lemma  3.1.  The  function 


m  y>  . 

f(Pji  Po»  •••»  P_)  =  f  6.  £n(l  -  e  )  +  (n.  -  6.)  p. 
14  m  j  1  J  J 

is  concave  for  0  <  p.  z  - 

Proof .  This  easily  follows  since 

| — I  =  -  6  •  ePl/(l  -  ePl)2  <.  0. 

V 

Since  f(pj,  •••»Pm)  is  concave,  and  the  constraints 
on  the  p^’s  are  linear,  the  following  lemma  follows  by 
arguments  similar  to  those  used  in  Theorem  1  of  Dykstra  and 
Madsen  (1974). 
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Lemma  3.2.  Assume  that  A  and  A  denote  column  vectors, 

£(•)  is  a  concave  real  valued  function  defined  on  an 
appropriate  subset  of  Rn ,  and  that  Ak  maximizes  f(A) 
subject  to  the  constraints  A|  A  <  b..  ,  i  =  1,2,  ...,k. 

Then  if  the  additional  constraint  Ak  +  j  A  <  is 

imposed,  we  may  assume  that 

a.  if  Ak+1  Ak  <  bR+1,  Ak+1  =  Ak; 

b.  if  Ak+1  Ak  >  bk+1,  A'+1  Ak+]  =  bk+1. 

The  key  in  determining  which  constraints  need  by 

imposed  in  maximizing  (3.2)  is  given  in  the  following  theorem. 

Theorem  3.1.  If  the  actual  restricted  M.L.H.'s  are 

expressed  in  the  form  of  (3.4) 

then  k,  >  k-  2  ...  2  k  *  2  0. 
i  l  m 

Proof.  Without  loss  of  generality,  assume  kj  <  k2>  Then 
if  the  maximum  is  found  with  the  m*  -  1  imposed  constraints 

X2  V 

1  Pj  =  C2’  **"  f  Pj  =  Cm*’ 

* 

the  value  of  k^  which  corresponds  to  the  first  constraint 

* 

is  such  that  kj^  <  <  ^2’  T^is  easily  follows  since  the 

Pj  are  nondecreasing  functions  of  the  kj .  Since 
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for  all  <  i  <  it  follows  that 
i  i 

I  p.*  ^  I  q.  for  i  *  1 . m, 

1  3  1  J 


where  the  p^.*  denotes 
imposing  only  the  m* 
this  implies  that  p  ^  * 
contradiction . 


the  values  which  maximize 
1  constraints.  By  Lemma 
=  Pj  which  is  an  obvious 


f(') 

3.2, 


Theorem  3.2.  The  algorithm  given  in  Section  II  obtains 

A 

Pj. 

Proof .  Straightforward  from  previous  considerations. 

The  closed  form  expression  for  a  particular  p. 
given  in  Section  II  can  also  be  obtained. 


Proof  (of  Theorem  2.1).  Clearly  the  expression  is  valid 
for  p^  by  Theorem  3.2.  Thus,  using  the  notation  of  the 
algorithm  given  in  Section  II,  it  easily  follows  that 

k,  =  t „  >  .  .  .  ,.  k  . 

1  1  -  2  -  m 


Let  i  denote  the  first  integer  <  i^  such  that  k. 


"1 

Z  q 

1 


ll  nj 

j  ’  f  tn(- 


6.  *  kj 


n.  +  k. 
J  1 


n.  -  6.  +  k. 
)  >  E  Ini-  3  3 


i-  1 
E 

1 


'“j  +  \ 


Then 


il  nj  ’  6j  +  ki 
I  £n( - ) 

*  "j  *  *, 


i-1 

2  E  q. 


which  is  a  contradiction.  Thus  no  such  i  exists  and 
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A 

=  k.  for  i  =  1,2,  .  ..,  i^.  Essentially  the  same 

arguments  handle  subsequent  intervals. 

We  should  note  that  whenever  <$j  =0,  (no  deaths 

are  observed  in  the  interval  [S^  ,  Sj  +  1))  ,  p.j*0,and  the  survival 

function  P(t)  places  zero  probability  over  this  interval. 

Heur  i  st  ica]  Iv,  the  stochast  i  ca  1 1  \  ordered  M.L.l:. 

acts  like  it  has  k.  additional  items  on  test  at  the  time 

J 

A  A 

S..  Over  the  interval  [S.  ,  S  .  ,  k.  -  k .  .  of  these  items 

J  J  J+1  J  J+1 

are  lost,  etc.  Note,  however,  that  the  need  not  be 
integers . 

IV.  THE  STOCHASTICALLY  ORDERED  TWO  SAMPLE  PROBLEM 

Let  us  now  consider  the  case  where  our  observations, 
perhaps  censored,  come  independently  from  two  different 
discrete  populations.  As  in  Section  III,  we  denote  the 
true  survival  functions  and  the  imposed  ordering  by 

P(t)  >  Q(t)  for  all  t. 

As  before,  we  assume  that  P(0)  =  Q(0)  =  1  WLOG . 

In  a  manner  similar  to  that  used  in  Section  III, 

the  function  to  be  maximized  may  be  expressed  as 

f(p,  q)  =  I  [6.  £n(l  -  )  +  (n.  -  6.)  p.]  + 

~  ~  ^  J  J  J  J 

[dj  fn (1  -  e^)  +  (m..  -  d..)  qi ] 

where  0  <  S,  <  ...  <  S  denotes  the  ordered  times  of 
1  m 

death  of  the  combined  sample,  (d..)  denotes  the  number 
of  deaths  at  S^.  in  the  first  (second)  sample,  n^  (nu) 
denotes  the  number  of  items  surviving  just  prior  to 
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in  the  first  (second)  sample,  and  p  (q^ )  represents  tn  f  P(S  j )  /P(  S  j_,)  I 

(£n  [Q(S  .  )/Q(S  .  _  ^)  )  .  The  order  constraints  are  still  of 
3  3 

the  form 

1  -i 

2  p .  >  2  q .  f  o  r  i  = 1 ,  ...»  m . 

1  3  1  3 

If  we  wish  to  maximize  f  subject  to  the  one 
constraint 

1  i 

2  p.  -  2  q.  , 

1  3  1  3 

i  i- 1 

then  we  may  let  q-  *  2  p.  -  2  q.,  and  set  the  partial 

1  1  3  1  3 

derivatives  equal  to  zero.  This  results  in  the  set  of 
equations 

p  -  p . 

-  6.  e3/(l  -  e3)  +  n. 

3  3 

•  q  • 

-  d.  eq  3  /  ( 1  -  e  3)  +  m. 

3  3 


-  di  e  7(1  -  e  *)  +m.  -  d.  ■  j  <  i . 


P.;  Pi 

-  e  3  (1  -  e  3)  +  nj  -  6^  = 

q . 

-  d.  eq3  (1  -  e  x)  +  m,  -  d.  = 

J  J  J 


-  5.  =  0 
3 


-  d.  =  0 
3 


J  '•  i 


eqi  /  (1  -  eqi)  -  m.  ♦  d±  =  k,  j 


<  i 
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Solving  these  equations  results  in  the  solutions 


n.  - 

6. 

A 

P3 

’  £n(-1HT 
3 

—4, 

3  > 

i 

m.  - 

d- 

A 

q3 

3 

-L), 

3  > 

i 

n .  - 

6.  + 

k 

A 

P3 

= 

3 

J.. 

+  k 

-), 

3  s  i 

m .  - 

d.  - 

k 

•n 

<cr 

* 

1 

-V 

3  s  i 

where  k  is  a  real  solution  to  the  equation 

i  n .  -  6 .  +  k  i  m.-d.-k 

?  3n.  +Jk  ■  >  -  ? 

1  J  1  3 

If  more  constraints  are  imposed,  similar  solutions 
can  be  obtained  for  disjoint  strings  of  the  and  q^ . 

Once  again  the  key  question  is  which  constraints 
need  be  imposed  in  finding  the  true  restricted  M.L.E.'s. 
However,  the  same  methods  and  lemmas  used  in  proving 
Theorems  2.1,  3.1  and  3.2  will  also  suffice  in  the  two- 
sample  case  and  lead  to  the  expressions  given  in  Section  TI. 

Another  appealing  method  of  handling  the  two-sample 
problem  which  will  work  if  the  samples  do  not  include  any 
censored  data  is  to  convert  it  to  a  one-sample  problem  as 
considered  in  Section  3  by  making  the  observation  that 
there  must  exist  a  survival  function,  say  R(t) ,  depending 
on  the  data  such  that 

P(t)  *  R(t)  a  Q (t)  for  all  t. 
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Intuitively,  the  more  equality  constraints  imposed 

A  A  A 

between  P(t)  and  Q(t),  the  more  P(t)  will  be  pulled  down 

A  A  A 

and  Q(t)  forced  up  (assuming  the  constraints  P(tJ  >  Q(t) 
are  already  imposed)  .  The  limiting  case  of  this  occurs 
if  all  equality  constraints  are  imposed  in  which  case 

P(t)  =  R(t)  *  Q  (t) 

A 

where  R(t)  denotes  the  Kaplan-Meier  (1958)  product  limit 
estimater  obtained  from  the  pooled  data.  Thus  if  we 

A  A 

construct  R(t) ,  treat  it  as  fixed,  and  compute  P(t)  as  in 

A  A 

Section  III,P(t)  will  be  our  restricted  M.L.E.  for  P(t)  . 

A 

With  obvious  changes,  Q(t)  can  be  obtained  similarly. 
Unfortunately,  this  method  does  not  work  when  the  samples 
include  censored  data  although  the  end  results  are  very 
close  to  the  actual  M.L.li.'s. 


In  the  special  case  of  no  censored  data  explicit 


expressions  for  the  k  .  are  given  by 

a  }  D 


b  b 

(I  <5^  ma  -  (E  dA)  na 
a  a 


b 

Zddi  *  5.) 
a  1 


This  is  of  course  a  weighted  average  of  hi  and  -n. 

iJ  ‘I 

with  weight  proportional  to  the  number  of  deaths  in  the 
two  populations  for  the  appropriate  interval.  In  this 
case  the  expressions  in  Section  2  are  equivalent  to  those 
in  Brunk,  et  al  (1966). 
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V.  AN  EXAMPLE 

To  illustrate  the  method,  we  apply  it  to  some 
data  gathered  by  Dr.  Martin  Alpert,  Department  of 
Cardiology  of  the  University  of  Missouri  Medical 
Center.  Dr.  Alpert 's  data  consists  of  survival  times 
for  people  who  have  had  heart  pacemakers  implanted. 

We  wish  to  estimate  the  survival  functions  separately 
for  males  and  females,  and  will  impose  the  constraint 
that  the  survival  function  for  females  never  drops 
below  that  of  males  since  it  is  well  documented  that 
females  are  longer  lived.  The  data  includes  many 
censored  observations  of  people  who  were  lost  to  the 
study.  The  data  on  pages  18  and  19  has  been  coded  to 
conform  to  the  notation  used  in  the  paper. 

The  M.L.E.'s  of  ordered  survival  functions  for 
males  and  females  is  shown  on  page  20.  KM-P(T)  and 
KM-Q(T)  denote  the  unrestricted  Kaplan-Meier  estimates 
for  the  females  and  males  respectively,  while  P(T)  and 
Q(T)  indicate  the  M.L.E.'s  obtained  using  our  algorithm. 
We  note  that  P(T)  has  been  forced  up  from  KM-P(T)  while 
Q(T)  has  been  forced  down  from  KM-Q(T) . 

To  try  and  get  an  idea  of  the  overall  effect  of 
our  order  restrictions,  we  computed  the  expected  values 
corresponding  to  our  various  survival  functions  when 


IWI  ~ 
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they  were  all  truncated  at  97  months.  These  expectations 
were 


P(T)  Q(T)  KM-P(T)  KM-Q(T) 

72.737  68.609  69.842  70.837  . 

Thus  we  see  that  our  order  restrictions  increased  the 
estimate  of  expected  life  (if  truncated  at  97  months)  by 
nearly  3  months  for  females  while  decreasing  it  by  approxi¬ 
mately  2  months  for  males 

A  computer  routine  (  in  Fortran)  for  implementing 
this  procedure  is  available  upon  request. 
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